Stability of localized modes in PT-symmetric nonlinear potentials 
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PACS 42 . 65 . Tg - Nonlinear guided waves 

PACS 42 . 65 . Sf - Optical instabilities (quantum optics) 

Abstract. - We report on detailed investigation of the stability of localized modes in the nonlin- 
ear Schrodinger equations with a nonlinear parity-time (alias PT) symmetric potential. We are 
particularly focusing on the case where the spatially-dependent nonlinearity is purely imaginary. 
We compute the Evans function of the linear operator determining the linear stability of local- 
ized modes. Results of the Evans function analysis predict that for sufficiently small dissipation 
localized modes become stable when the propagation constant exceeds certain threshold value. 
This is the case for periodic and tanh-shaped complex potentials where the modes having widths 
comparable with or smaller than the characteristic width of the complex potential are stable, while 
broad modes are unstable. In contrast, in complex potentials that change linearly with transverse 
coordinate all modes are stable, what suggests that the relation between width of the modes and 
spatial size of the complex potential define the stability in the general case. These results were 
confirmed using the direct propagation of the solutions for the mentioned examples. 
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Introduction. — Since the introduction of the con- 
' cept of parity-time- (VT-) symmetric potentials [II, non- 
I Hermitian Hamiltonians possessing purely real spectrum 
' have received considerable attention [5] due to their rele- 
vance to the quantum mechanics and optics. In the con- 
text of optical applications, it was natural that the concept 
. was generalized to the nonlinear systems ^ , where the ex- 
' istence of localized modes was shown to be possible (we 
notice that experimental observation of the T'T'symmetry 
in linear optics was recently reported [4]). As a natural 
extension of this activity, the existence of stable localized 
nonlinear modes in nonlinear PT-symmetric lattices has 
been recently demonstrated in [5]. Nonlinear gain and 
loss compensating each other were also addressed recently 
within the framework of the nonlinear dimer model 6 . 
A general important property of nonlinear T'T-symmetric 
systems is that they admit continuous families of local- 
ized modes parameterized by the propagation constant, 
similarly to what happens for the nonlinear Hamiltonian 
systems. This is in spite of the need to satisfy the balance 
between dissipation and gain, as this happens in dissipa- 
tive systems of general type where localized modes appear 
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as attractors, rather than elements of a family of solu- 
tions. In this sense, the models with T^T-symmetric po- 
tentials occupy a special "place" between the Hamiltonian 
and dissipative systems. 

A striking effect related to the existence and stability of 
the localized modes in nonlinear T'T-lattices is that they 
can be stable even in the absence of modulation of con- 
servative part of nonlinearity [5]. The stability, however 
depends on the relation between the mode width and the 
period of the potential. More specifically only sufficiently 
narrow modes were found to be stable. Now we observe, 
that in the limit when the width of a mode goes to zero 
and the potential is a smooth function, the behavior of 
the imaginary part of the potential can be approximated 
by the respective linear function of the coordinate (recall, 
that the imaginary part is an odd function). For instance, 
the purely imaginary nonlinear PT-potential i sin(277) con- 
sidered in [5] can be approximated as i sin(2r7) ~ 2iri for 
sufficiently narrow modes strongly localized about rj = 0. 
This leads to the natural question: is it possible to obtain 
stable localized modes in a nonlinear T^T-symmetric po- 
tential of a general type, provided the widths of the modes 
are small enough? Here we give an affirmative answer to 
this question and illustrate stable modes for three quali- 
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tatively different case examples; of a periodically varying 
dissipation and gain, i.e. the lattice case, of dissipation 
and gain tending to constants at the infinity (see ([6]) be- 
low), and to the linearly increasing gain and losses (the 
model (ESI). 

The aim of the present paper is the analytic study of 
the linear stability and dynamics of solitons in nonlinear 
T^T-symmetric potentials, giving positive answer to the 
above question. The main model we are interested in is 
the complex nonlinear Schrodinger equation [5]: 
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^d'^q- [I + V ill) +iW{i^)]\q\'q 



(1) 



In the optical applications q is the dimensionless elec- 
tric field propagating along the ^-direction (^ > 0) with 
?7 (ry g R) being the transverse coordinate. We no- 
tice that the physical mechanism for both gain and loss 
are well known. The former one refers to the standard 
two-photon absorption, which becomes dominating mech- 
anism of losses in semiconductors at sufficiently high in- 
tensities, while nonlinear amplification can be realized in 
electrically-pumped semiconductor optical amplifiers (see 
e.g. [7]). We are interested in the special case of a VT- 
symmetric nonlinear potential where the V{rii) and W{ff) 
are both real and obeying the relations: 



V{ri) = V{-r]) and W{ri) = -~W{-Ti). 



(2) 



Localized modes. — We look for stationary local- 
ized solutions of eqs. ([I])-©, which can be searched in 
the form q(^, rj) = w{r])e'^^^ subject to the boundary con- 
ditions lim^_j.±oo |9(Cj ^)l = 0- Bearing in mind optical 
applications, we refer to b as to the propagation constant. 
The stationary wave function w{'if) obeys the equation 



-d^w - 6w + [1 + y (r;) + iWi-q)] \w\^w = 0. 



(3) 



Let us calculate how many parameters one has to intro- 
duce in order to unambiguously define a localized mode 
w{r]). To this end, let us first agree that we identify 
the modes w{ri) and 'w{ri)e^^ , € R, which are not dis- 
tinguishable from physical point of view. Then a spe- 
cific symmetry of eq. (jS]) induced by relations ([2]) sug- 
gests that without loss of generality the localized mode 
w{ri) = Writ]) -\-iWi{ri) can be chosen to be T^T-symmetric, 
i.e. having even real part and odd imaginary one: 



Wr{r]) = Wri-r]), Wi{r/) = -~Wi{-r]). 



(4) 



Let wiji) be some solution of eq. (jS]) vanishing as ry — > 
-|-oo, but not necessarily vanishing as 77 — >■ —00. Then 
for 77 — > +00 the nonlinear term in eq. ([3]) is negligi- 
ble and in the corresponding limit w{ri) is described by 
the linear equation ^d^w — bw — 0. Thus for 77 +00 

the solution w(r/) behaves as w{7j) = Ce*'^[e"^'' -I- o(l)], 
where C and ip are real constants. For a generic solu- 
tion 'w{ri) the constant (p is not zero, but let us tem- 
porarily restrict ourselves to the case ip = 0. Eqs. (|1]) 



dictate that if the solution w^rj) represents a localized 
mode, then 10(77) must obey d,jWr{0) = and 7«i(0) — 0. 
Now, let us admit nonzero ip in the asymptotics for w{r]). 
Obviously, this leads just to multiplication of ?Z'(77) by 
the factor e^'^. Therefore, we can formulate a weaker 
condition for 711(77) to represent a localized mode: there 
must exist such ip that Wi{0) cos 95 -1- 70^(0) sin 1^9 = and 
driW r (0) cos (fi — 9^77)^ (0) siu iy9 = 0. These equations are 
compatible if and only if 



(5) 



Thus any localized mode 77;(77) can be identified with a so- 
lution of the eq. ([S]) which contains two real unknowns: C 
and b. If we fix one of them, typically this is the propaga- 
tion constant b, then eq. ([5]) results in one or several solu- 
tions for the parameter C what indicates that eq. ([3]) ad- 
mits continuous families of localized modes for fixed V{'r]) 
and 1^(77). This feature is typical for 7^7~-symmetric (lin- 
ear or nonlinear) potentials and constitutes significant dif- 
ference compared to the conventional dissipative systems. 

In fig. [1] (a) we show two families of localized modes 
plane {b,U), where U = J \q\^dri is the energy flow (the 
integration limits are omitted wherever the integration is 
over whole real axis) , obtained for the nonlinear potential 



^"(77) = 0, and W{t]) ^ crtanh(27/) 



(6) 



with a — 0.7. The modes exist only if b exceeds certain 
threshold value. The modes of the lower (upper) curve 
in fig. [1] (a) can be referred to as fundamental (higher) 
modes. A typical profile of a fundamental mode is shown 
in fig. [1] (b). Below we focus on the fundamental modes. 

On the other hand, the localized modes can be consid- 
ered as bifurcating from the limit cr = 0, where eq. ([1]) re- 
duces to the conventional nonlinear Schrodinger equation. 
In fig. 1 (c) we show a branch of fundamental modes on the 
plane (cr, U) found for a fixed value of b. The continuation 
from the limit ct = will be used below as an approach 
for analytical investigation of stability of the fundamental 
modes. 

Linear stability analysis. — Substituting the per- 



turbed solution q(^, 7/) 

where \p±\ <C |77;|, and the overline stands for the com- 
plex conjugation, in eq. ^ and linearizing it around 77;(7;) 
one arrives at the eigenvalue problem Lp = Xp, where 
p = (75+ +p-,i{p- —73+))"^ (hereafter the superscript T 
stands for matrix transposition) and L is given by 



L 



192 

2^V 



N21 



2 V ' 



(7) 



where 



Nil = -2[1 + V{r])]wrW^ - W{r])[3wl + wf], (8a) 

N22 - 2[1 + Virj)]wrw, - W{rj)[wl + 3w^], (8b) 

Ni2 = b-[1 + V{7^)][wf + 3w^] - 2W{r])wrW,, (8c) 

N21 = -b+[l + y(77)][37«2 + w^] - 2W{T^)WrW,. (8d) 
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Fig. 1: Properties of the modes for the potential ([6]). (a) 
Families of the fundamental and higher locahzed modes for 
a — 0.7. (b) Real and imaginary parts Wr^i as well as the mod- 
ulus u = \w\ of the fundamental mode at 6 = 0.35, a — 0.7. (c) 
The branch of the fundamental modes found for b = 0.35. The 
circle corresponds to the mode shown in panel (b). (d) Real 
part Xr = ReA of perturbation growth rate vs a for b = 0.2 
(curve 1) and 6 = 0.35 (curve 2). 



The mode w{ri) is unstable if and only if there exists an 
eigenvalue A with positive real part. 

For V{ri) — W{ri) = the localized mode w{ri) is a 

standard NLS soliton: Writ]) — wf'^'' — \/26sech(\/26?7) 
and 'Wi{r]) — wf^'' = 0. Designating the operator L in 
this case by L'-^\ we recall that the spectrum of L*^"^ is 
well known [5]. In particular, the point spectrum of L'-'^' 
consists of the only eigenvalue Aq = which is isolated 
and has algebraic multiplicity (a.m.) equal to 4 and geo- 
metric multiplicity (g.m.) equal to 2. The eigenfunctions 
corresponding to Aq read 

V'^°^^(a,«.(°),o)^ ^i°) = (o,^(°))^ (9) 

There also exist two generalized eigenfunctions, namely 



0, 



and ■022'' 



(10) 



such that Lip^^j 



I (0) ■ 



1,2. Here dbwi'^'' is obtained 
„(") 



by means of differentiation of Wr with respect to b. 

Let us now assume that V{r]) and W{ri) are not equal 
to zero but satisfy T'T-symmetry relations ([2]). If at the 
same time V{ri) and W{i]) are small enough, then they 
can be considered as a perturbation to the operator L'"^. 
Behavior of the spectrum of L'"' subject to the pertur- 
bation determines linear stability of the localized mode. 
In particular, it is crucial to understand what happens to 
the multiple eigenvalue A = since a generic perturbation 
of the operator L^^^ leads to splitting of the eigenvalue 
A = into several simple eigenvalues. As a result, unsta- 
ble eigenvalues can arise in the vicinity of A = 0. 




(b) 
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Fig. 2: Possible behavior of the Evans function E{X) for real A 
is schematically illustrated, (a) d\^E > and apart from the 
zero Ao = the Evans function has two real zeros ±Ai. (b) 
dl^^E <0 and no real zeros of the Evans function exists besides 
of Ao = 0. Instead, there are two purely imaginary zeros ±Ai. 



The multiplicity of the eigenvalue A = is related to 
rotational (i.e. phase) and translational symmetries of 
the model. The dissipative perturbation introduced by 
V{j]) and W{ri) breaks the translational symmetry and 
preserves the rotational one. Due to the last fact A = 
remains an eigenvalue for L. However ip^^ ceases to be 
the eigenfunction corresponding to A = and the only 
eigenfunction for A = is given by ■012 = {—Wi,Wr)'^. 
Disregarding parity of the functions Wr^i{ri) we observe 
that if A is an eigenvalue of the operator L then A is also 
an eigenvalue. Then, recalling eqs. (|4]) we find that —A 
is an eigenvalue, as well. Therefore upon the perturba- 
tion exactly two simple eigenvalues arise in the vicinity of 
A = and those eigenvalues have opposite signs and either 
purely real or purely imaginary (see fig. 

Evans function for T^T-symmetric potentials. — 

The operator L can be associated with the Evans function 
E{X) [HIlO], which is an analytic function defined on the 
whole complex plane except for the points of the continu- 
ous spectrum of L. An important property of the Evans 
function is that it has a zero at some A if and only if A is 
an eigenvalue of L. In addition, the order of that zero is 
equal to the a.m. of the eigenvalue A. 

Before proceeding with explicit definition of the Evans 
function for the case at hand, we recall that perturba- 
tion of the operator L leads to that the eigenvalue Aq has 
a.m. = 2. Hence the Evans function corresponding to the 
perturbed operator L has a zero of the second order at 
A = 0: E{0) ^ dxoE ^ 0, dl^Ey^O (hereafter d^^ stays 
for jth partial derivative with respect to A evaluated at 
A = 0). Without loss of generality one can assume that 
E{X) < for all A > 1. Then the stability of the station- 
ary mode is determined by the sign of d^^E: if df^E > 
then E{X) necessarily has exactly one positive and one 
negative zero what corresponds to instability; vice versa, 
if 9? -E < then both roots of the Evans function lie on 
the imaginary axis, and the solution is stable. These con- 
siderations are illustrated in fig. [H 

In order to define the Evans function for the operator 
L given by eq. (O we rewrite the eigenvalue problem [L — 
X)p = in the form of the system of four first-order ODEs 
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driY — MY, where Y = {pr,Pi, d^Pr, driPif^ and 



for aU a the fohowing equahties hold: 



M 



( 


-2iV2i 
V 2(iVn - A) 






2(A - 
2Ni2 



1 0\ 

1 



0/ 



(11) 



The coefheients Njk are given by eqs. ([S]). Parity of the 
functions Wrirj) and Wi{r]) imposed by eqs. (|H) implies that 
N12 and N21 are even while Nn and N22 are odd as func- 
tions of 77. 

For the sake of simplicity we suppose that W{t]) — 
aWQ{ri), where cr ^ 1 is the small parameter while 
V{r]) = 0{a^). Then the coefficients Njk become also 
dependent on a. Obviously, now substitution a — )■ — a is 
equivalent to substitution rj ^ —rj what in turn is equiva- 
lent to Wr(r]) — )■ Wrirj) and Wi{r]) — —Wi{r]). As a result, 
N12 and A^2i are even while Nn and N22 are odd in a. 

Now we choose four solutions Yj{ri; a; X), j — 1,...4, 
such that Yi and Y2 are linearly independent in rj and 
vanish as 77 — >■ —00 while I3 and I4 are linearly in- 
dependent and vanish as 77 — > +00. Then the Evans 
function is given as 4 x 4 77-independent determinant [3] : 
E{a;X) = det[Yi,Y3,Y2,Y4]. In order for E{a;X) to be 
unambiguously defined and to depend analytically on a 
and A, we explicitly fix the choice of the solutions Yj set- 
ting 



^i,3(^;0;0) = 
1^2,4(r/;0;0) = 



(0,u;(°),0,9,«;(°))' 



(recall that rp^^ and ■0^2^ are given by eqs. ([9])). Requiring 
this choice to be consistent with parity of the coefficients 
Njk we find that 

Y3i7T,±a;TX)^ Y,{-7]; a; \) * J± (12a) 
Y4(r/;±a,TA) =±l2(-?7;a;A)* J± (12b) 

where J+ = (-1, 1, 1,-1)^, J_ = (-1,-1, 1, 1)^, 
and * stays for the element- wise multiplication of the ma- 
trices. Eqs. (|12l) imply that the Evans function is even 
with respect to both its arguments: E{(t; A) — E{a; —A) — 
E{—a; A). Respectively, the second derivative of the Evans 
function at A = Aq can be searched in the form of the fol- 
lowing expansion: 



diE{a;\o)^a''E2k + oia"'), E2k ¥^ 0, 



(13) 



where /c = 1, 2, ... is to be determined. 

Subject to the perturbation, i.e. for a 0, the eigen- 
function ipn disappears; the solutions Yi^3(77; cr; 0) become 
unbounded and no longer correspond to any eigenfunction 
of the operator L. However, the eigenfunction 1^12 defined 
above persists and depends smoothly on cr. As a result. 



Y2Airi;T;0) 



■012 

dn4'i2 



{-Wi, Wr, -dr,Wi, dr,Wr) 



Since Aq has a.m. =2 in the spectrum of the perturbed 
operator. Lemma 3.3 from [10] can be applied. It states 
that d\Y2{'q; a; 0) = 9a 14(77; a; 0) for all 77 and a. Bearing 
in mind these facts and differentiating straightforwardly 
the Evans function one arrives at the following expression 
for the second derivative of the Evans function at A = 0: 



dlEin;^) - A^X\Yr,Y3,d\X^2 ~ Y^),Y2\. 



(14) 



Let us now recall that E{a; A) is an even function of a. 
It means that da^d\^E = where d^^ stands for the jth 
partial derivative with respect to a evaluated at (Tq = 0. 
Calculating the second derivative with respect to a and 
evaluating it at <tq one arrives at 

dldlE = Aet[dlXYi - I3), 1^1, 9^,(1^2 - I4), I2] + 
det[a,„(yi - Y3),Y^,d,,dl^{Y2 - 1-4), I2]. (15) 

Using parity of the coefficients Nij and the symmetries of 
the solutions Yj given by eqs. one can recognize that 
r.h.s. of eq. (jlSp generically is not equal to zero, i.e. we 
can set fc = 1 in eqs. (fT3| and obtain: 



(16) 



d„dlE{Q;Q)=0, dldlE{Q;Q)^Q. 



We failed to obtain expressions for the mixed deriva- 
tives d1^{Yi — Y3) and daQd\^{Y2 — I4) which would al- 
low for efficient analytical or numerical computation of 
d^^df E{0;0). However, using the information given by 
(1161) one can use another representation for derivatives 
of the Evans function [10] which fits better for analytical 
and numerical investigation. To this end, we recall that 
apart from the eigenfunction ■0i2 there exists a generalized 
eigenfunction 022 such that L1IJ22 = ■012- The adjoint op- 
erator for L reads = L^. There also exists the adjoint 
eigenfunction v such that L^v = 0. It is found in [10] that 
for any cr 7^ the second derivative of the Evans function 
at A = can be found as 9^£'(cr;0) = 2{xp2i,v), where 
(•, •) represents standard inner product of vector- valued 
functions. Taking into account eqs. (|16p one can construct 
an expansion for ('02i, v) with respect to the small param- 
eter cr, < |cr| ^ 1. To this end, we firstly write down an 
expansion for the stationary mode itself. Eq. ([3]) dictates 
that this expansion acquires the following specific form: 



(0) 



+ a^wi^Uo{a^), 



aw. 



(1) 



+ o{a'^), where the 



coefficients wl?^ and 



„(i) 



solve the equations 



L+U;(2) = ) - Wo{v)wi% (17a) 

L-w\''> =~Wo{v){wi'^f, (17b) 

(2) 

Notice that wi- ' is an 



andi± = p2 



d'^-b+i2±l)(wi''^y 



even function of?/ while w^ is odd. Next, using eqs. (jlOp 



p-4 



Localized modes in PT-symmetric nonlinear potentials 



and the definition of ■022 (i-e. the equation = 1P12) 

we obtain an expansion for the generalized eigenfunction: 



■022 = 









.(1) 



+ o(a), (18) 



where the coefficient x^^"^ is an odd function of r/ solving 
the equation 

= -wi"'>dbwi"'>[3Wo{v)w^°^ + 2wf ^] + ^ (19) 

Using the definition the adjoint eigenfunction v and re- 
quiring the derivative d\E{(j] 0) to satisfy constrains ([1 
we observe that expansion for v must have the form: 



V = a 







drjWr 



(0) 



•,(1) 



Then p^^^ is given as p 



(1) 



f3w. 



(0) 



(20) 

+ /, where / is a par- 



ticular solution of the equation 

L-f = w^°^d„wi"'> [2u;f ' - Wo{v)wi% (21) 

The solvability condition for eq. ((2T|) (i.e. orthogonality 
of its r.h.s. to wl^^) is automatically provided as long 
as eq. (|17bp holds. At the same time, the coefficient /3 
should be chosen to satisfy the solvability condition of the 
equation with respect to 0^^): 

d^wi°^6wlP'>wi^^ + (wf - 2Wo{v)wi"^w^% (22) 

Respectively, we require r.h.s. of eq. (1^^ to be orthogonal 
to dr^wi^^ what yields (3 = I1/I2, where 



and h = — J djjWo{ri) (w^.^^^^ drj. Finally, using eqs. ([T 

and (1201), we can rewrite eq. ((T^ in the form: d\E{a\ 0) = 
2{ip2i, v) — (7^ E2 + o{a'^) , where the coefficient E2 is given 

asE2^2 j {dr,w'^^\^^^ +af,u;(°V^^^) dri. Functions x^^^ 

and p^^-* can be computed numerically from the linear 
equations and (pij) . what gives an algorithm for ob- 
taining the coefficient E2 ■ For given b positive value of E2 
corresponds to the situation when the modes are unstable 
for small a. Vice versa, negative E2 implies stability of 
the modes for small a. 

Discussion of the results and Conclusion. Let 

us now turn to the results of the stability analysis of the 
particular examples (see fig. [3]). We start by recalling the 
results for the case Vijf) = and W{rf) = asm{2rf) re- 
ported earlier in [5]. Physically, this case corresponds to 
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0.0 
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Fig. 3: Domains of existence and stability of fundamental lo- 
calized modes on the plane {b,a) sin-shaped (a), tanh-shaped 
(b), and linear (c) potentials W{rf) and V{Tf) — 0. (d): The 
dependence of the coefficient E2 on b in sin-shaped (curve 1), 
tanh-shaped (curve 2), and linear (curve 3) potentials. 



the pure dissipative nonlinear lattice where domains with 
the nonlinear gain alternate with the nonlinear dissipa- 
tion. This case was already discussed in [5], and in par- 
ticular it was shown that for sufficiently small a the non- 
linear modes become stable if the propagation constant 
b exceeds a threshold value ^^''(ct). The developed here 
approach based on the analysis of the Evans function al- 
lows us to compute numerically 6'^'"(0), as this is shown 
in fig. 13] (d) (the black curve 1) from which one observes 
that the coefficient E2 changes its sign at b^'^{0) « 1.05 
which corroborates results reported in the panel (a). This 
critical value corresponds to narrow modes, whose widths 
are of order of the period of the dissipative lattice. 

Next we consider potential (jB]) [see also fig. [1] which is 
characterized by only one domain with nonlinear gain and 
one domain with nonlinear dissipation. From fig. |3] (b) 
we observe that the domains of existence and stability are 
similar to those obtained for the sin-shaped potential: in 
particular, there exists the threshold for the stability of 
the modes. Quantitatively, now the stable modes can be 
broader than in the case of periodic potential: threshold 
value of propagation constant is remarkably lower than in 
the previous case. More specifically, using the approach 
based on the Evans function we obtain [the blue curve 
in panel (d)]: 6^''(0) ~ 0.49. This result agrees with 
fig. |3] (b) and with fig. [1] (d) where dependence of per- 
turbation growth rate ReA on a is shown. 

From the results for sin- and tanh- shaped potentials we 
can conjecture that (i) the most stable modes are localized 
on the scale where the dissipative potential can be approx- 
imated by the linear function, i.e. where W{r]) ~ 2ar] on 
the width of the mode; and (ii) properly introduced gain 



p-5 



D. A. Zezyulin et al. 




Fig. 4: Propagation of initially perturbed nonlinear modes of 
eq. (HSl) for a = 0.5 > a""" (a) and <j = 0.2 < a"" (b). For both 
the cases 6=1. 

and dissipation from the both sides of the mode enhances 
its stability. These arguments readily lead to the model 

idiq^-\dlq-\qW-2iaT^\qW (23) 

where all modes in the limit ct — >■ should be stable. This 
is indeed what happens, as one can see from fig. [3] (c): the 
instability threshold 6^''(0) disappears completely. This is 
also confirmed by analysis based on the Evans function 
[the red curve 3 in the panel (d)]. 

Now from panel (c) we observe that both lines, i.e. up- 
per border of existence domain cr^^^ (black curve) and the 
upper border of stability domain a'^^ (red curve), follow 
the parabolic scaling law b ^ . This law can be un- 
derstood from the simple scaling arguments as follows. If 
q{^,'q) is a solution of eq. (^5)) . then g(^,?7) = ^^(Ci^) 
with ^ — fj^f and fj — arj is a. solution of the complex NLS 
equation without any parameter: id^q — —^d^q — \q\'^q — 
2ifi\q\^q. This observation as well as the fact that a nar- 
row mode "feels" only the local dissipative term, allows 
one to make further conclusions about the behavior of the 
curves in fig. |21 Since the mode widths tend to zero at 
6 — >• oo (and a bounded), and all the examples of W{r]) 
were chosen to have the same slope at 77 = 0: W{r]) ~ 2ari, 
the existence curves a^^^ in the panels (a) and (b) tend at 
5 — > 00 to the parabola shown in the panel (c). This 
conjecture was supported by numerical simulations up to 
b — 20 where for the linear, tanh-, and sin-shape poten- 
tials we found (t"pp — 3.18, 3.20, and 3.19, respectively. 

In fig. [4] we present two examples of evolution of unsta- 
ble (panel a) and stable (panel b) modes. These results, 
obtained by direct integration of eq. (|23l) [i.e. of the par- 
ticular case of the model ([T|) which is mostly "exposed" 
to eventual nonlinear instabilities] confirm the results of 
linear stability analysis. Thus, the modes belonging to 
the stability domains propagate undistorted over indefi- 
nitely long distances, even if they are strongly perturbed 
initially. 

To conclude, we have investigated fundamental modes 
in imaginary T'T'-symmetric nonlinear potentials. Such 



potentials allow for existence of localized modes which 
are stable at least in the limit when the mode is narrow 
enough. The stability was established both as the linear 
stability, on the basis of the Evans-function analysis, and 
using direct numerical study of the mode evolution (notice 
that while the direct propagation ensures also the nonlin- 
ear stability of the modes in a finite domain, the nonlinear 
nature of the perturbation may introduce new features of 
the nonlinear stability of the solutions on the whole real 
axis). Although our analysis was performed for nonlinear 
potentials, the established symmetry properties have more 
general character and the approach can be applied also for 
linear 'PT-potentials, as well as to the cases where both 
linear and nonlinear 'P7~-symmetric potentials are present. 
In the latter case stability of the modes may change dra- 
matically. For example, in presence of a periodic linear 
PT-potential broad small-amplitude modes are expected 
to be stable as long as imaginary part of the linear poten- 
tial is below a critical value. This situation is in contrast 
to the one shown in fig. [3] (a) where broad modes are un- 
stable. As another interesting question, we would like to 
mention the exploration of asymmetric nonlinear modes 
similar to ones reported in [B], although their existence 
may require more sophisticated nonlinearity landscapes. 

^ ^ ^ 
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